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ABSTRACT 


We  discuss  a  fundamental  radar  imaging  problem  for  narrowband  electromagnetic  waves  that 
extends  the  recent  results  originally  obtained  in  [1,  2]  in  the  scalar,  or  acoustic,  setting.  In  particular  we 
demonstrate  the  ability  to  efficiently  image  three-dimensional  convex  conducting  bodies  by  using  the 
knowledge  of  the  scattered  electric  field  for  one  fixed  monochromatic  illumination  of  the  target.  In  this 
problem  our  knowledge  of  the  scattered  electric  wave  is  understood  to  be  in  the  form  of  measurements  of 
the  amplitude  and  phase  of  the  tangential  components  of  the  radiated  electric  field  (generated  by  the  radar 
target)  on  a  finite  two-dimensional  array.  We  also  provide  a  more  general  version  of  this  streamlined  result 
that  describes  imaging  of  multiply  connected  nonconvex  objects  with  the  same  measurements. 
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1.  INTRODUCTION 


We  begin  with  a  few  introductory  remarks  that  outline  what  we  mean  by  narrowband 
electromagnetic  imaging.  The  problem  concerns  estimating  the  shape  and  size  of  an  object  on  the  basis  of 
known  information  in  the  form  of  measurements  of  reflected  waves  taken  at  some  distance  from  the 
unknown  object.  Specifically,  how  do  we  form  a  three-dimensional  image  of  a  potentially  very  distant 
object  using  something  like  a  radar  or  collection  of  such  sensors?  Figure  1  depicts  such  a  measurement 
scenario. 


Figure  1.  Data  collection  scenario. 


The  problem  is  simple:  image  the  scatterer,  or  target,  with  a  single  incident  wave  with  the  aid  of 
measurements  of  the  scattered  wave  taken  on  a  finite  two-dimensional  array.  For  example,  the  scattered 
wave  is  a  frequency-dependent,  i.e.,  narrowband,  induced  voltage  response  created  by  the  illumination  of 
the  target.  Although  the  underlying  physics  and  mathematics  that  lead  us  to  the  imaging  algorithm  can 
appear  complicated,  the  algorithm  we  will  ultimately  employ  is  extremely  easy  to  understand  and 
implement. 

The  imaging  algorithm  may  be  summarized  as  follows.  The  measured  electric  field,  as  a  function  of 
two  variables  on  the  array  and  for  a  chosen  origin  of  the  problem,  is  decomposed  into  its  generalized 
Fourier  modes,  and  a  combined  one-dimensional  power  spectrum  of  this  signal  is  generated.  The  critical 


1 


point  in  the  power  spectrum  where  the  combined  spectral  energy  transitions  from  a  highly  oscillatory 
behavior  and  begins  a  super  exponential  decay  is  then  found.  This  index  is  determined  by  computing  the  /2 
norm1  of  the  accumulated  energy  of  these  modes  and  defining  the  cutoff  at  ninety  percent  of  the  total 
spectral  energy.  The  index  cutoff  is  converted  to  a  length  through  division  by  the  wave  number  (271/^)  of 
the  narrowband  waveform  used  in  probing  the  target.  This  length  represents  the  radius  of  the  smallest 
sphere,  centered  at  the  origin  we  previously  specified,  which  must  contain  the  unknown  scatterer.  The 
process  is  then  repeated  with  new  center,  or  origin,  and  terminated  when  enough  centers  of  interest  have 
been  considered.  Finally,  the  intersection  of  all  the  computed  spheres  is  computed  and  the  resulting 
common  region  is  taken  as  an  approximation  to  the  convex  scatterer. 

We  mention  that  the  power  spectrum  possesses  three  very  pronounced  and  particular  behavior 
regimes  as  a  function  of  increasing  mode  number.  Figure  2  illustrates  the  three  regions  of  interest  and 
shows  the  combined  power  spectra,  labeled  as  |  (3n  |  for  a  candidate  example.  The  first  regimes  takes  place 
within  0  <  n  <  45,  the  second  roughly  within  45  <  n  <  55,  and  the  final  region  for  all  n  beyond  this. 


1. 


i.e.,  L2 


norm  of  a  discrete  collection  of  complex  numbers. 
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Figure  2.  Oscillation  and  decay  of  the  combined  power  spectra. 


First,  there  is  a  region  of  highly  oscillatory  modes,  followed  by  a  region  of  extremely  rapid  decay, 
finally  followed  by  a  region  of  essentially  null  values.  The  index  at  which  this  power  spectrum  begins  this 
rapid  decline  defines  the  radius  of  the  smallest  sphere  with  a  specified  origin,  which  must  contain  the 
scattering  object. :  By  considering  several  origins  we  construct  several  such  spheres,  each  of  which  must 
contain  the  scatterer.  Hence,  by  considering  the  intersection  of  all  such  spherical  regions  we  arrive  at  an 
approximation  to  a  three-dimensional  image  of  the  assumed  target.  Figure  3  provides  an  illustration  of  this 
process. 


2.  For  the  moment,  we  mention  that  the  reason  for  this  is  intimately  related  to  the  unique  analytic  continuation  of  the 
measured  electric  field  on  an  array.  We  will  revisit  this  in  the  section  treating  the  supporting  theory  behind  our  main 
result  and  defer  further  discussion  until  such  time,  due  to  its  technical  nature. 
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Figure  3.  An  example  of  two-dimensional  projection  of  the  intersection  of  three  spherical-containing  regions  and  the 
true  convex  target. 


The  remainder  of  this  report  reviews  the  physics  that  govern  this  problem  and  how  that  physics  and 
some  mathematical  analysis  lead  to  a  well-understood  and  fully  justified  target  characterization  algorithm. 
Moreover,  this  treatment  suggests  additional  innovative  object  characterization  techniques  through  its 
description  of  the  fundamental  scattering  processes  at  work.  These  results  will  be  presented  in  forthcoming 
articles. 
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2.  THE  PHYSICAL  MODEL 


We  begin  with  a  brief  overview  of  the  fundamental  inverse  scattering  problem  at  the  heart  of  our 
application.  Our  imaging  technique  will  employ  novel  processing  of  measurements  of  electromagnetic 
waves.  Hence,  to  illustrate  the  origins,  success,  and  limitations  of  the  algorithm,  we  begin  with  a  brief 
review  and  characterization  of  what  we  term  the  forward,  or  direct,  scattering  problem:  given  a  well- 
defined  scattering  object  and  a  specified  incident  wave,  we  formulate  the  response  to  that  incident 
illumination  and  analytically  predict  how  such  a  wave  would  interact  with  the  scatterer. 

We  start  with  the  time-dependent  Maxwell’s  equations  in  free  space,  i.e.,  in  a  subdomain  of  Q  c  R 
in  which  there  are  no  electric  or  magnetic  current  densities,  for  which  the  electric  field  £  and  the  magnetic 
field  B  satisfy 


V  a  £(x,  t)  +  p03(B(;t,  r)  =  0, 


(1) 


V  a  B(x,  t)  —  e03,8(.r,  t)  =  0. 


(2) 


Here,  x  =  (xx,  x2,  *3)  e  /?3  is  our  spatial  variable  and  (e  R  denotes  time.  The  constant  quantities 
p0  and  e0  are,  respectively,  the  magnetic  permeability  and  electric  permittivity  of  free  space.  We  briefly 
mention  that  the  electric  and  magnetic  fields  are  complex-valued  vector  quantities  as  functions  of  both 
position  and  time,  so  that  £,  B:  R3  x  R— »C3.  Additionally,  V  a  is  the  standard  curl  operation  and  d , 
denotes  differentiation  in  time. 

Rather  than  dealing  directly  with  the  problem  in  the  time  domain,  we  will  recast  this  system  of  six 
coupled  partial  differential  equations  into  one  that  depends  on  frequency.  On  the  basis  of  this  system,  we 
arrive  at  our  narrowband,  or  fixed-frequency,  imaging  result. 

To  accomplish  this  goal,  we  employ  the  Fourier  transform  mapping  functions,  or  distributions,  in 
time  into  ones  depending  on  angular  frequency  w.  In  this  way  we  define  frequency-dependent  electric  and 
magnetic  fields  as 


E(x,w)  =  je  '"'Z(x,t)dt  and  H{x,w)  -  Je  B(x,  t)dt , 
R  R 

whose  inversions  are  given  by 


(3) 
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(4) 


E(x,t)  =  —  \e  E(x,w)dw  and  B{x,  t)  =  —  \e  H(x,w)dw. 
m J  2n J 

R  R 

We  note  that  i  =  4- 1  represents  the  imaginary  unit. 

Employing  the  Fourier  transform  essentially  converts  differentiation  in  time  into  multiplication  in 
frequency.  Hence,  we  arrive  at  the  frequency-dependent  Maxwell  system 


V  a  E(x,  w)  -  iw\i0H(x,  w)  =  0,  (5) 

V  a  H(x,  w )  +  iw£0E(x,  w)  =  0  .  (6) 

We  now  take  a  moment  to  distinguish  between  two  fundamentally  different  types  of  electric  and 
magnetic  fields.  By  this  we  mean  outgoing,  or  scattered  or  radiating  waves  (E5,  Hs),  and  incoming,  or  local 

incident  plane  waves  (£',  Hl).  We  define  the  total  fields  we  would  measure  with  any  physical  device  as  the 
sum  of  these  two  types  of  waves.  Namely,  we  will  use,  and  so  we  define  here,  the  total  fields  E  and  H  as 


E(x)  =  E\x)  +  E?{ x),  H(x)  =  H\x)  +  H\x) .  (7) 

Any  outgoing  pair  of  electric  and  magnetic  fields  must  satisfy  the  so-called  Silver-Miiller  radiation 
condition,  which  we  recall  is  given  by  either 


lim  (Hs  Ax-rEs)  =  0,  r  =  |jc| 


(8) 


or 


lim  ( E?  ax  +  rtf)  =  0,  r  =  |jc| 


(9) 
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Equivalently,  Es  =  (E\,  E?2,  Ef3)  and  Hs  =  (H\,  H2,  H\)  means  that  each  Cartesian  coordinate 

v  3  5  3 

component  Ey.R  — >  C  and  Hy.R  =  1,2,3,  of  the  two  vector  fields  E?  and  Hs  satisfies  the 

Sommerfeld  radiation  condition,  which  is 


lim  ( dr-ik)E  =  lim  ( drik)H  =  0,  r  =  |x|  .  (10) 

r  — >  oo  J  r  — >  co  J 


The  Silver-Muller  and  Sommerfeld  radiation  conditions  imply  the  asymptotic  behaviors 


«>W 

E\x)  =  -j-^-(£'00(i)  +  0(1*1"  )),  (11) 

and 

ik\x\ 

H\x)  =  e—{HM  +  0(\4 -‘)),  W->oo,  (12) 

where  x  =  (x/|x| )  and  |x|“  =  x]  +  x2  +  x\ .  Here,  the  quantities  E^  and  denote  the  far-field  patterns, 
or  scattering  amplitudes,  of  the  outgoing  waves  Es(x)  and  Hs(x).  We  emphasize  that  for  most  remote 
electromagnetic  sensing  experiments,  these  are  the  quantities  we  measure.  Additionally,  the  electric  and 
magnetic  scattering  amplitudes  defined  on  the  unit  sphere  with  normal  v  2  must  also  satisfy  the  asymptotic 
boundary  conditions 


=  <  V,  HJ  =  0,  where  //„  =  VfJ  a  E„.  (13) 

We  mention  that  this  equation  simply  means  that  the  electric  and  magnetic  scattering  amplitudes  are 
purely  tangential  vector  fields,  possessing  no  normal  contribution  in  any  direction  x  e  S2  . 

Incoming  electric  and  magnetic  waves  are  ones  similar  to  the  radiating  case,  only  they  satisfy  the 
conjugate  of  the  Silver-Muller  and  Sommerfeld  radiation  conditions.  That  is,  rather  than  having  constant 
phase  fronts  that  migrate  outward,  in  the  usual  sense,  they  travel  and  converge  on  the  origin  as  a  function 
of  positive  time.  Additionally,  by  a  local  plane  wave  propagating  in  the  direction  y  we  mean  a  certain  type 
of  incoming  wave — specifically,  a  Herglotz  wave  with  a  Dirac  kernel  supported  in  the  direction  y . 

We  will  now  turn  our  attention  to  the  case  of  the  scattering  of  an  incident  wave  pair  ( E H1)  from  a 
perfect  conducting  scatterer.  In  particular,  what  we  mean  is  that  the  tangential  components  of  the  total 
electric  field  must  vanish  throughout  the  surface  of  the  scatterer,  so  that  if  vD(x)  denotes  this  normal  of  the 
scatterer  D  at  some  point  x,  then  we  have 


7 


vD(x)  a  E(x)  =  0,  Vjcg  dD. 


(14) 


In  summary,  the  problem  is  to  understand  how  an  incident  pair  of  fields  (E,  H1)  interacts  with  the 
conductor  D  and  generates  the  scattered  pair  of  fields  ( Es,  Hs)  as  a  function  of  various  parametrizations  of 
the  incident  wave  energy,  e.g.,  their  frequency  or  energy  and  direction.  We  note  that  the  problem  may  be 
stated  and  treated  as  a  system  of  partial  differential  equations,  admitting  very  insightful  analysis  and 
solution  techniques  via  the  theory  of  finite  elements.  For  more  on  this  approach  see  [3, 4],  For  the  purposes 
of  the  exposition  of  this  report,  however,  we  will  choose  to  frame  this  problem  in  terms  of  integral 
equations  for  the  clarity  they  provide  in  terms  of  the  operator  mappings  between  the  scattering  system 
inputs  (E‘,  H‘)  and  its  outputs  (Es,  Hs). 

In  any  domain  Cl  c  R  within  which  there  are  no  charge  or  current  sources,  each  type  of  outgoing  or 
incoming  wave  satisfies  our  frequency-dependent  Maxwell  system 


Va£i,!-  ikH,  s  =  Va  H''s  +  ikEl's  =  0,  k  =  w/c  =  w(|i0£0)  ',  (15) 


and  hence,  the  total  fields  E  and  H  do  as  well.  This  system  of  coupled  partial  differential  equations  has  an 
alternative  integral  equation  representation  called  the  Stratton-Chu  integral  equation,  which  for  a  perfect 
conductor  has  the  form 


E(x)  =  E\x)  +  a  V  a  $  G(x-y)v(y)  a  H(y)ds(y) ,  (16) 

dD 

H(x)  =  H\x)  +  V  a  jG(x-y)v(y)  a  H(y)ds(y).  (17) 

dD 

Here,  ds(y)  denotes  surface  measure  on  the  boundary  of  the  conducting  scatterer  D,  and  the  integrals 
exist  as  improper  ones  as  x  — >  dD  3  y .  Again,  we  emphasize  that  the  pair  ( E ,  H)  denotes  the  total 
electromagnetic  and  magnetic  fields,  and  we  identify  the  integral  kernel  G(x-y)  appearing  in  both  integral 
equations  as  the  fundamental  solution  of  the  scalar  Helmholtz  equation.  Specifically,  G  is  the  outgoing 
fundamental,  or  distributional,  solution  of 


3 

(A  +  k2)G(x-y)  =  8(x-y),  x,ye  R 3,  A  =  X  (18) 

j=  i 
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where  8  e  £'(/?  )  is  the  so-called  Dirac  distribution,  or  “delta  function,”  and  is  given  by 


G(x-y) 


1  g,fcU-yl 

4 it  U-yl  ’ 


x^y . 


(19) 


The  asymptotic  behavior  of  the  convolution  kernel  G  is  of  extreme  importance;  hence,  we  take  a  few 
moments  to  examine  it.  Suppose  the  point  y  is  contained  within  the  finite  ball  of  radius  R  centered  at  the 
origin,  which  we  write  as  BR(0).  Suppose  further  that  x  =  x/\x\  g  S2  =  {x  g  /?3:  x]  +  x\  +  x\  =  1 }  and 
G(x)  is  the  outgoing  fundamental  solution  to  the  operator  A  +  k2  in  R3.  Then, 


G(x-y) 


1  g ikM 

4n  |jc| 


-«'*  <■*.  y) 


x  — >  oo  . 


(20) 


This  claim  is  evident  when  we  examine  the  asymptotics  of  the  term  |.t  -  y|.  First,  if  y  is  required  to 
remain  in  a  finite  region,  while  x  may  tend  to  infinity,  then 


-y\  =  W  +  +  o(\x\  2), 


(21) 


Since  y  e  ^^(0) , then 


max 

£g  f\(x,y)\  </?<-•  (22) 


Hence, 


i  ik\x\ 

G(x-y)=  ^~-e~,k{x,y)(l+0(\4~')),  W-»~.  (23) 

The  asymptotic  behavior  of  the  acoustic  outgoing  fundamental  solution  G  coupled  with  the 
asymptotic  behavior  of  the  scattered  electric  and  magnetic  waves  implies  that  the  electric  and  magnetic 
scattering  amplitudes  are  given  by 


£«,(•*)  =  j  e~'k{x,y)  v((y)  a  H(y))  Axds(y) 

ao 


(24) 
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(25) 


Hco(x)  =  ~XA  je  lk(x,y)v((y)  AE(y))ds(y). 

dD 


The  formulations  for  the  electric  and  magnetic  scattering  amplitudes  now  provide  us  an  operator  we 
term  F,  whose  analytical  properties  will  be  of  significant  interest.  It  tells  how  and  why  we  can  image 
conducting  bodies  with  narrowband  measurements  of  the  scattered  electromagnetic  waves  made  at  a 
distance.  For  this  reason  we  will  examine  a  few  more  aspects  of  the  operator  F  and  briefly  discuss  how  it 
indicates  and  describes  the  interaction  between  an  incident  incoming  magnetic  wave  H1  and  a  conducting 
body.  More  importantly,  we  will  hit  on  how  such  a  description  benefits  us  in  suggesting  a  reconstruction, 
or  three-dimensional  imaging  algorithm. 

From  the  previous  analysis  we  see  that  the  electric  scattering  amplitude  is  generated  by  the  surficial 
current  j  =  vdaH  generated  on  the  boundary  of  the  conductor  D  by  the  interaction  of  an  incident  pair  of 
waves  ( E H‘) .  In  summary. 


Ecc(x)  =  (Fj)(x):  =  J  K(x,y)j(y)ds(y),  (26) 

dD 

where  the  kernel  of  operator  F  is 


K{x,yY  =  -^e-‘kCx'y)XAXA.  (27) 

471 

If  we  recall  and  write  the  second  of  the  Stratton-Chu  integral  equations  as  a  function  of  the  incident 
magnetic  wave 


H'a  =  jeik{y,x)a (x)ds(x),  a  eTL2(S2),  (28) 

S2 

where  the  subscript  indicates  the  parametrization  of  H ' — e.g.,  in  this  instance  H'a  is  a  superposition  of 
entire  plane  waves,  also  called  Herglotz  waves,  and  plays  the  role  of  the  weighting  of  such  square 
integrable  tangential  vector  waves  and  is  the  Herglotz  kernel — as 


H(x)  =  H'a(x)  +  V  a  jG(x  — y)v(y)  a  H(y)ds(y) 

dD 


(29) 
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=  :  H'(x)  +  (BH)(x), 


(30) 


then,  formally  given  an  known  parametrization  of  an  incident  magnetic  wave  H'a ,  the  electric  scattering 
amplitude  is  generated  by 


EM)  =  J  E(x,y)vD  a  ((I  -  B)~l  H‘a)(y)ds(y)  (31) 

ao 


=  :(F  H‘a)(x),  ieA. 


(32) 


In  the  next  section,  we  present  the  main  result  of  this  report,  which  is  that  the  coefficients  in  the 
generalized  Fourier  series  representation  of  £„  on  A  encode  the  boundary  of  the  scatterer  dD.  Moreover, 
this  result  gives  us  a  robust  and  efficient  algorithm  that  can  reconstruct  the  scatterer  with  limited  but 
appropriately  sampled  values  of  E „  on  A.  The  details  of  the  properties  of  F  and  the  representation  of  the 
signal  E ^  on  the  array  A  that  support  these  results  will  then  be  discussed  in  detail  in  Section  4. 
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3.  MAIN  RESULTS 


Before  discussing  the  main  results  we  need  to  address  two  issues.  The  first  concerns  the  dependence 
of  the  electric  scattering  amplitude  on  the  origin  of  the  physical  problem.  Since  we  are  at  liberty  to  define 
the  origin  where  we  choose  and  this  value  changes  the  values  of  the  measured  quantity  we  have,  we  must 
determine  this  influence.  Moreover,  we  will  demonstrate  that  this  analysis  leads  us  to  a  very  attractive  and 
easily  implementable  algorithm  for  characterizing  the  scattering  body.  The  second  issue  regards  the  vector 
representation  of  the  tangential  vector  field  on  our  measurement  array. 

We  begin  with  a  quick  treatment  of  the  continuous  dependence  of  the  electric  scattering  amplitude 
with  the  origin  c  e  R '  .  In  what  we  developed  in  the  previous  section,  we  assumed  the  origin  to  be  the 
standard  one,  namely,  c  =  0.  However,  if  instead  we  choose  any  arbitrary  point  c  to  be  the  center  of  our 
universe,  then  the  electric  scattering  amplitude  is  generated  by 


EM)  =  ( Fj(-c))Cx ) 


(33) 


=  J  E(x,  y)j(y  -c)ds(y)  (34) 

dD-c 

=  “4I  1  e~‘k{X'y)x^xA  j(y-c)ds(y)  (35) 

dD  -  c 


ik 

471 


ik(x,  c ) 

e 


J  e  ,k{x’y)x  ax  Aj(y)ds(y) 


dD 


(36) 


=  O  cCx)(Fj)(x), 


Oc.(x)  =  e 


ik(x ,  c > 


(37) 


where  dD-c  is  the  translate  of  dD  to  the  ocentered  origin.  Hence,  given  the  true  scattering  amplitude 
produced  by  the  physical  interaction  with  the  incident  wave,  we  may  produce  its  value  if  instead  we  insist 
that  the  point  c  is  the  origin  by  performing  a  unitary  operation,  i.e.,  a  simple  phase  multiplication  at  each 
point  where  we  measure  the  field.  We  further  add  that  this  concept  of  shifting  and  choosing  various  origins 
plays  a  very  powerful  and  useful  role  in  the  imaging  algorithm. 

We  recall  that  the  electric  scattering  amplitude  we  measure  in  remote  sensing  on  a  flat  or  curved  two- 
dimensional  array  will  consist  of  two  orthogonal  components.  By  this  we  mean  that  the  electric  scattering 
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amplitude  is  a  vector-valued  function.  Also,  since  any  flat  or  smoothly  curved  array  is  locally  equivalent  to 
the  unit  sphere,  we  begin  with  a  very  brief  discussion  of  two-dimensional  vector  fields  supported  on  the 
unit  sphere  S2.  In  what  follows,  for  a  scalar  function  i/(0,  <p)  defined  on  the  unit  sphere  S2,  we  will  use  the 
so-called  surficial  gradient,  defined  by 


1  du  du\ 
sin  03(j>’  dds  ’ 


(38) 


which  is  an  element  of  the  tangent  space  of  the  unit  sphere,  denoted  by  . 

Let  v  be  a  complex-valued  vector  field  supported  on  the  unit  sphere  S2.  Then  we  say  that 
v  e  TL  (S  ) ,  provided  that  in  the  representation  of  the  vector  field  on  the  sphere  S2 


=  I  I 


xeS 


n  =  1  m 


(39) 


obeys  the  condition 


|(|<v(jc),  Vs2l/"(i))  +  (v(x),xa  V^Y^ix))  ^ds(i) 


<  oo 


(40) 


where  the  terms  (v,  and  (v,  x  a  represent  the  projections  of  the  tangential  field  on  S2 

onto  the  orthonormal  bases  V^2 Y"  and  xaV^Y"  .  Here,  the  functions  Y"(x)  are  the  spherical 
harmonics,  of  order  n  and  degree  m  on  5 2  with  x  =  sin0,  cost)),  sin0sin<|>,cos0,  0  e  [0, 7t]  and 
<)>  €  [0,  2jc]  .  For  our  purposes  we  define  the  harmonic  functions  in  the  following  way: 


Oe,  <i>) :  = 


Myiu-lmhy,,  os  ,, 

V  4;t  (n  +  |m|)! 


„  2n  +  1  (n  -  m)\  n,n,  n.  i 

em  \— \ - , - COS0)e 

m*'  4n  (n  +  m)!  n 


-lm,  m  >  0 
1,  m  <  0 


(41) 
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Above,  the  functions  P™( cos0)  are  the  associated  Legendre  functions;  c.f.  [5,  6]  for  more  details  on  their 
various  properties  and  governing  equation.  This  definition  of  the  spherical  harmonics  allows  us  to  express 
their  complex  conjugates  as  simply 


Cm- 


(42) 


Also,  this  definition  yields  a  complete  orthonormal  set  on  L~(S a  in  the  sense  that 


J KCx)K\i)dSW  -  C"'.  (43) 

s2 

where  8”*’ is  the  Kronecker  delta.  Lastly,  we  note  that  the  finiteness  criterion  amounts  to  requiring  that 
the  generalized  Fourier  coefficients  satisfy 


z  z  kf +|p:|2<oo-  <44) 

n  =  1  m  =  -n 

We  now  have  all  the  key  ingredients  in  place  and  state  the  main  result  of  this  article  as  the  following 
theorem. 

3 

Theorem  1 .  Let  BR(c)  denote  the  sphere  of  finite  radius  R  centered  at  the  point  c  e  R ‘  ,  and  let  E ^  be 
known  on  some  section  T  cz  of  the  unit  sphere  and  have  a  series  expansion  in  tangential  spherical 
harmonics  defined  on  T, 


ESx)  =  Z  X  +  ier,  (45) 

n  =  1  m  =  -n 

with  Fourier  coefficients  depending  on  R  and  c.  Then  the  following  are  equivalent 

1 .  Br(c)  contains  dD 

2.  There  exists  a  magnetic  current  j  supported  within  BR(c),  depending  on  R  and  c,  such  that 


£«(*)  =  J  K(x-y)j(y)dy ,  xg  T. 

Br(c) 


(46) 
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3.  The  weighted  energy  of  the  Fourier  coefficients,  depending  on  R  and  c,  is  finite. 


oo  n 

I  "„(*.<»  I  |oC(K,<f +  |p;(R,c)|2,  (47) 

n  =  1  m  -  -n 

Remark  1.  The  theorem  states  several  important  facts  all  at  once.  In  summary,  each  item  implies  one 
another,  which  is  to  say  that  given  any  item,  the  other  two  are  true  as  well.  In  what  follows  we  address 
these  implications  and  provide  some  explanation. 

To  begin  it  states  that  if  the  scatterer  dD  can  be  placed  or  fully  contained  within  some  sphere  of  radius  R 
centered  at  the  point  c,  then  there  exists  a  magnetic  current  supported  on  the  boundary  of  a  two- 
dimensional  surface  contained  within  this  sphere,  which  when  integrated  against  the  kernel  of  the  electric 
scattering  amplitude  generator  yields  the  known  electric  scattering  amplitude  £„ . 

Next,  the  result  states  that  such  a  scattering  amplitude  must  have  partial  Fourier  coefficients  that  decay 
faster  than  (2 n/ ek{R  +  |c|))  ” ,  as  n  tends  to  infinity.  This  is  the  observable  quantity  we  exploit  for  various 
centers  c,  which  then  determines  the  value  of  R  for  each  such  c.  Upon  intersecting  all  such  spheres  for  an 
appropriate  collection  of  centers,  this  quantity  yields  an  approximation  of  the  convex  hull  of  the  scatterer, 
which  if  the  scatterer  is  convex  is  an  approximation  of  the  scatterer  itself.  We  note  that  in  the  limit  that  |c| 
tends  to  infinity  and  the  approximations  of  the  convex  hull  of  dD  become  better  and  better,  theoretically 
speaking.  From  a  numerical  point  of  view  this  approximation  becomes  more  and  more  difficult,  as  it 
requires  the  computation  of  higher  and  higher  Fourier  modes,  which  is  difficult  due  to  the  rapid  oscillation 
of  the  spherical  harmonics  for  large  n.  However,  even  numerically  speaking,  this  approximation  can  still  be 
done  if  we  utilize  arbitrary  machine  precision  techniques. 

Lastly,  if  what  we  call  the  weighted  energy  of  the  Fourier  coefficients  converges  for  some  R  and  center  c, 
then  it  must  be  that  there  exists  a  sphere  with  such  a  radius  and  center  that  contains  the  totality  of  the 
scatterer  dD,  which  in  turn  implies  that  the  electric  scattering  amplitude  is  generated  by  some  magnetic 
current  on  the  surface  of  some  two-dimensional  manifold  contained  within  such  a  sphere. 

In  a  sense,  what  is  going  on  here  has  to  do  with  the  analytic  continuation  of  an  analytic  function.  The 
function  £ro  is  analytic  in  several  variables,  x  and  k,  for  instance.  Furthermore,  it  is  analytic  in  x,  which 
means  if  we  know  it  on  some  distant  surface  away  from  the  scatterer  dD,  the  function  may  be  analytically 
continued  up  to  the  boundary  of  the  scatterer.  Where  this  continuation  is  no  longer  possible  coincides  with 
R  for  any  given  c,  and  means  that  no  longer  will  the  weighted  Fourier  coefficients  converge,  and  no  longer 
is  there  a  magnetic  current  that  can  be  integrated  over  such  a  smaller  sphere  and  give  rise  to  the  known 
field  £„ . 

We  now  offer  a  concise  version  of  this  result  in  the  form  of  an  algorithm  that  has  been  successfully 
implemented  in  the  acoustic  setting. 

Algorithm  1.  (Convex  3D  Imaging). 

1 .  Measure  E ^  on  some  two-dimensional  array  A. 
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2.  Determine  the  tangential  analysis  eigenfunctions  spanning  the  array  A  that  are  used  to  determine 
the  spectrum,  i.e.,  generalized  Fourier  coefficients,  of  the  signal  E ^  (functions  such  as  V  and 
x  A  VC  in  the  case  when  A  is  the  unit  sphere3). 

3.  Optimally  span  a  collection  of  points  c  and  compute  the  associated  generalized  Fourier 
coefficients  for  that  center. 

4.  Determine  the  radius  R  for  each  such  center  by  examining  where  the  partial  sums  of  the  Fourier 
coefficients  begin  a  rapid  transition  to  zero. 

5.  Intersect  all  such  spheres  to  define  the  approximation  of  the  convex  hull  of  the  scatterer. 

Remark  2.  We  are  also  compelled  to  make  another  remark  concerning  the  uniqueness  of  reconstructing  the 
true  scatterer  in  this  problem.  Our  main  result  states  that  the  measured  asymptotic  wave  we  know  may  be 
analytically  continued  in  a  unique  fashion  up  to  some  natural  boundary  that  corresponds  to  the  boundary  of 
some  scatterer.  To  be  fair,  we  must  mention  that  in  general  this  scatterer  need  not  coincide  with  the  true 
scatterer  of  interest.  Rather,  what  we  ultimately  arrive  at  in  this  process  is  the  so-called  scattering  support 
of  the  measured  field,  or  the  scatterer.  See  [1, 2]  for  more  of  the  fundamental  details  concerning  this  object, 
as  this  object  can  be  highly  nontrivial  and  rather  complex.  For  our  purposes  here,  we  define  this  object 
through  the  following  convention. 

Definition  1 .  Let  ch(supp  f)  denote  the  convex  hull  of  the  support  of  a  distribution / and  let  s  =  -1/2.  Then 
the  convex  scattering  support  of  the  far  field  E ^  for  a  perfectly  conducting  electromagnetic  scatterer  is 


cSksupp  Ex  -  p  ch(suppf), 
Ff  =  E_ 


(48) 


where  /e  TL2(/?(/?3))  4 


What  we  can  uniquely  determine  is  the  common  element  to  all  candidate  scattering  objects  that 
could  give  rise  to  the  measured  field  of  interest.  We  need  to  make  this  interjection  and  state  this  result, 
since  to  date  there  is  no  rigorous  proof  that  a  perfectly  conducting  scatterer  may  be  uniquely  determined 


3.  For  more  complex  measurement  arrays,  these  eigenfunctions  may  be  analytically  or  numerically  computed  by 
treating  a  standard  eigenvalue,  eigenfunction  problem  whose  formulation  follows  from  the  spectral  theory  of  compact 
linear  operators  presented  in  the  next  section. 

4.  The  space  of  distributions  /^(Z?3)  for  5  =  -1/2  is  the  standard  Sobolev  space  jhT]/2(R3)  of  distributions  having 
compact  support  in  i.e.,  they  are  TL 2  vector  fields  on  a  two-dimensional  manifold  such  as  the  surface  of  a  sphere. 
Moreover,  since  we  are  assuming  that  the  scatterer  is  perfectly  conducting,  we  may  specify  s  =  -1/2  rather  than  having 
to  take  more  general  values  of  5  falling  within  the  interval  (-2,  0)  for  which  the  original  theory  of  scattering  support 
was  developed. 
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from  one  single  fixed  energy  excitation  with  measurements  taken  on  such  a  section  of  the  unit  sphere,  or 
two-dimensional  array.  Nonetheless,  on  the  basis  of  uniqueness  results  and  several  numerical  studies  of 
this  problem  in  the  acoustic  setting,  it  is  widely  suspected  that  the  convex  scattering  support  and  the 
scatterer’s  true  support  are  one  and  the  same  in  the  case  of  a  perfectly  conducting  electromagnetic  obstacle. 
A  careful  study  and  proof  of  this  conjecture  would  constitute  a  major  breakthrough  in  the  academic  realm 
of  scattering  theory;  however,  it  is  beyond  the  current  scope  of  this  report. 

Remark  3.  We  should  also  take  a  moment  to  offer  a  generalization  of  this  result,  which  allows  us  to  image 
more  complex  three-dimensional  targets  that  need  not  be  convex.  In  the  case  of  multiply  connected 
nonconvex  targets,  we  can  offer  a  theoretically  based  method  that  can  image  such  systems.  However,  a 
practical  numerical  implementation  of  such  a  result  has  yet  to  be  done  and  is  beyond  the  work  presented  in 
[7]. 

More  generally,  we  have  the  following  formulation. 

Theorem  2.  If  Q.  denotes  a  (possibly)  multiply  connected  finite  subset  of  R3,  and  Em  is  known  on  some 
section  T  c  S~  of  the  unit  sphere  and  has  a  series  expansion  in  tangential  spherical  harmonics  defined  on 

r, 

oo  n 

£.(i)=S  S  cCvCw icr.  ,49) 

n  =  1  m--n 

with  Fourier  coefficients  depending  on  Q  then,  the  following  are  equivalent: 

1 .  Q.  contains  dD. 

2.  There  exists  a  magnetic  current  j9  supported  within  and  depending  on  Q,  such  that 


E~(x)  =  |^(i-y);(y)^y,  xe  r.  (50) 

a 

3.  The  weighted  energy  of  the  Fourier  coefficients,  now  depending  on  Q,  is  finite. 


o°  n 

2>.<n>  S  k(Q)|2  +  |u:(n)f.  (5i) 

n-  1  m  =  -n 

For  the  purpose  of  brevity  we  shall  not  go  into  the  rather  complex  nature  of  the  exact  asymptotics  of 
these  more  generalized  weights  wn(Q),  since  their  description  is  so  highly  connected  with  the  geometry 
and  number  of  components  of  Q.  We  shall  have  to  content  ourselves  that  this  formulation  extends  the 
previous  theorem  and  represents  the  more  general  methodology  of  imaging  nonconvex  scatterers,  since  the 
test  region  Q  is  allowed  to  be  nonconvex,  whereas  the  aforementioned  spheres  were  not. 
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In  the  supporting  theory  to  follow  in  the  next  section  we  give  the  appropriate  indications  that 
corroborate  the  claim  made  here  and  tie  it  to  a  result  previously  obtained  in  the  scalar  acoustic  scattering 
case  [7].  Additionally,  we  provide  the  mathematical  machinery  that  yields  a  numerically  feasible  algorithm 
that  can  image  such  scatterers  of  interest. 
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4.  SUPPORTING  THEORY 


We  begin  with  a  basic  physical  treatment  of  the  scattering  problem  that  inspires  the  more  rigorous 
supporting  theory,  which  underpins  the  two  theorems  given  in  the  previous  section. 

The  scattering  of  a  fixed  monochromatic  electromagnetic  wave  of  wave  number  k  =  2n/X  from  a 
perfect  conductor  dD  situated  in  the  ball  BR(0)  is  equivalent  to  the  radiation  problem 


V  a  Es(x )  -  ikH\x)  =  V  a  H'\x)  +  ikE'\x )  =  0,  x  e  R?\(dD  n  D) , 


(52) 


where 


E'Ud  =  ve  TL\dD) 


(53) 


is  specified.  The  general  solution  of  such  a  radiation  problem  in  the  complement  of  BR(0)  admits  the  form 


E\x)  =  X  Z  a>  A  ^\kr)Ynn{x))  +  p:v  A  V  a  (xh^\kr)V^x)),  r  =  W  >  R,  (54) 


v"'/ 


n  =  1  m  =  -n 


where  h(n' }  is  the  spherical  Hankel  function  of  the  first  kind — c.f.  [6] — and  is  defined  by 


h{n'\z)  =  jn(z)  +  iyn(z)  ,  (55) 

with  jn  and  yn  respectively  denoting  the  standard  spherical  Bessel  functions  of  the  first  and  second  kind. 

Examining  the  asymptotic  behavior  of  the  outgoing  electric  field  leads  to  the  electric  scattering 
amplitude,  which  is  known  to  have  the  form  [8] 


£_<;>  =  1  I  r'-’^v^i-oCvOA 


1  -  n 


(56) 
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Moreover,  the  tangential  electric  scattering  amplitude  is  an  element  of  TL2(S2)  only  if 


s(^)2"  I  l<|2+lfCr<”'  VK>diam(dD).  (57) 

n  =  1  tn  -  -n 

The  necessity  of  the  weighted  summability  of  the  Fourier  coefficients  a and  (3”'  is  a  consequence 
of  the  asymptotic  behavior  of  the  spherical  Hankel  functions  as  a  function  of  its  index  n.  For  a  fixed 
positive  argument  z. 


n 


n  — >  °°  . 


(58) 


Hence,  has  a  finite  L 2  norm  on  the  surface  of  any  sphere  of  radius  r  >  R  in  TL2(Br(0 ))  only  if  the  partial 

sum  of  the  combined  Fourier  coefficients  decay  faster  than  ’(Ul  )|2  as  n  ->  °°  for  all  r  >  |z| .  For  each 

such  r,  where  we  see  this  convergent  behavior  indicates  where  the  radiation  problem  may  be  posed  as  a 
homogeneous  Maxwell  system,  i.e.,  how  far  in  toward  the  origin  the  analytic  scattered  field  may  be 
uniquely  continued  as  a  free  solution  to  the  Maxwell  system.  That  value  r*  where  the  analytic  continuation 
may  no  longer  be  performed  defines  the  boundary  of  the  scatterer.  The  analysis  to  follow  treats  the 
operators  involved  in  the  scattering  problem  and  formalizes  this  basic  physical  principle. 


(59) 


To  begin  this  analysis,  we  note  that  we  need  know  only  the  electric  scattering  amplitude  on  some 
sector,  or  neighborhood,  of  the  unit  sphere.  The  significance  of  this  is  that  there  exists  a  collection  of 
coordinate  charts  { <py } ,  :/?“—>  5  that  cover  the  sector  of  interest  S2.  What  this  means  is  that  there  is  a 
change  of  coordinates  so  that  a  “piece”  of  the  unit  sphere  may  be  made  to  look  like  a  finite  two- 
dimensional  plane,  and  vice  versa.  Hence,  in  talking  about  the  measurement  array  A  or  the  unit  sphere  S2 
we  are  talking  about  the  same  object,  up  to  a  nonlinear  change  of  variables,  so  that  all  the  norms  to  follow 
can  be  made  equivalent  on  either  the  more  easily  treated  unit  sphere  S  or  the  true  measurement  array  A. 
See  chapter  two  of  [4]  for  more  information  concerning  this  differential  geometry. 

We  proceed  with  a  characterization  of  the  operator  that  maps  induced  magnetic  currents  on  the 
boundary  of  the  scatterer  to  the  electric  scattering  amplitude. 
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Proposition  1.  Let  dD  c  R2  be  smooth  surface  such  that  T(W  is  defined  everywhere,  and  suppose 
there  are  no  nontrivial  solutions  to  the  interior  version  of  the  Stratton-Chu  integral  equations.  Then 
F  :  TL'(dD )  — >  FL“(S“)  is  compact  and  linear  and  has  dense  range  in  TL2{S2) . 

Proof.  Linearity  is  obvious  as  F(g,  +g2)  =  Fgx  +Fg2  for  any  g,  2in  TL2(dD).  F  is  compact,  since  it  has 
a  so-called  Hilbert-Schmidt  kernel;  i.e.,  K  e  TL2 (ds(dD)  x  ds(S2)) .  Let  ge  TL2(dD)  be  a  magnetic 
current  on  dD.  Let  6  and  tp  denote  the  tangential  unit  vectors  on  S2.  Then,  the  previous  claim  follows, 

since 


\\F\\ 


TL2(ds(dD)xds(S2)) 


=  J  J  j9D^(*>3')^0')  2^W, 


(60) 


a 


- ik(x , 


1 6lt  S2  dD 


’fdsCx). 


(61) 


<-^pOD)p(52)<oo,  (62) 

16n‘ 

where  p(dD)  and  p(S2)  are  the  finite  surface  measures  of  the  scatterer  and  the  unit  sphere,  respectively. 
Now,  given  that  K  has  a  Hilbert-Schmidt  kernel,  we  find  that  given  g  e  TL2(dD)  a  magnetic  current  on 
dD, 


2  2  2 
II  ^gll  -  II  *11  TL2(ds(bD)  x  ds(S2)) ^  TL2(dD)  <  °°  ’ 

2  2 

proves  the  image  of  g  is  included  in  TL  (S“) . 

Lastly,  we  note  that  our  operator  F  is  equivalent  up  to  a  tc/2  rotation  of  the  tangential  field  Ex  on  S 2 

with  the  far-field  operator  presented  in  [8]  in  theorem  7.4,  which  is  known  to  have  dense  range  in 

2  2  .  7 

TL~(S  ) ,  provided  that  k  is  not  a  Dirichlet  eigenvalue  for  the  interior  problem  on  dD. 

We  now  address  the  range  of  the  operator  F  acting  on  magnetic  currents  on  the  surface  of  the  sphere 
Br(c)  in  terms  of  the  restriction  of  the  original  operator  F  to  TL2  current  distributions  supported  on  the 
boundary  BR(c).  Below,  F\dB^  denotes  this  restriction. 

Proposition  2.  Let  F|a^(c.)  be  the  restriction  of  F  to  BR(c).  Then  possesses  a  Hilbert 

adjoint,  such  that 
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F\1b«)  :  FLV)  ->  TL\dBR{c)) 


(64) 


and  for  g  e  TL2(S2) 

(Fl dB'{c)8)(y)  =  Avfi#(c),  yedBR(c).  (65) 

Proof.  Since  F  is  compact  and  linear  between  two  Hilbert  spaces,  then  by  the  Riesz  representation  theorem 
there  exists  a  bounded  linear  operator  F  coinciding  with  the  action  of  the  adjoint  F*  of  F  such  that  it  maps 
TL  ( S  )  to  TL  (dBR(c)) .  F  may  be  computed  by  equating  the  pairings 

S')  j-^2(S2)  ~  d*  F  TL2(dBK(c ))  (66) 

for  any  j  e  TL2(dBR(c))  and  g  e  TL~(S2) .  Hence,  we  find  that 

p(j\  8)  =  <Fj,g) rL2(52)  (67) 

=  f  <f  e~lk{x,y)xAXA j(y)ds(y),g(x))  ids(x)  (68) 

V  J3B*(c) 

=  J  <f  (e~'k{x,y)x  ax  a  j(y),  g(xj)ds(x))  ids(y)  (69) 

=  ^(y),-{vBR(c)^\^e~'kUy)i  ax  a  g(i)ds(x)j  a  vBr{c))  ^ds(y)  (70) 

=  0>  P  S)TL2(dBR(c))-  (71) 


24 


Hence,  the  claim  is  proved. 

Next,  we  address  the  so-called  singular  system,  or  singular  value  decomposition  of  F  restricted  to 
TL 2  currents  on  the  boundary  of  the  sphere  BR(c).  Since  F  is  linear  and  compact  then  we  know  it  admits  the 
general  representation 


oo  n  2 

*W)=Z  I  ©o>.O9>,c)0^.c).  (72) 

n  =  1  m  =  -n  j  =  1 

Here,  the  superscripts  (j  =  1 )  and  (j  =  2)  refer  to  the  tangential  components  in  the  9  and  <)>  directions, 
respectively,  and  the  term  ®  is  the  tensor  or  outer  product  of  the  two  orthonormal  bases  \\i^ \n  and  <p©m  for 
j  =  1,2.  We  add  that  the  action  of  acting  on  some  TL2  current  j  is 

oo  n  2 

(FlaW)(i)=  Z  Z  ®  °n!m(R’C)<J’Vnlm(R’C))Vn!m-  (73) 

n  =  1  m  -  -n  J  ~  ^ 

We  also  mention  that  this  decomposition  then  goes  to  define  the  Hilbert  adjoint  of  F|*a^(C.)  through 


oo  n  2 

F\*dB„(c)  =  Z  Z  ©  (74) 

n  =  1  m  =  -n  J  =  1 

where  the  above  overline  denotes  complex  conjugation.  Furthermore,  we  add  that  the  countable  singular 
system 


r  ^,0)  A>(j)  A,  O')-,00 

{<*»..«  .  m  >{„=!,  m  =  -„. 


(75) 


2  2 

of  F:  TL~(dBR(c))  — >  TL  ( A )  any  finite  two-dimensional  measurement  array  A  may  be  constructed  by 
determining  those  functions  such  that  any  of  the  following  identities  are  satisfied  for 


7=1,2. 


FI  „A0)  _  A,U).A,(j) 

I  dBK(c)^n,  m  ® n,m  V n,m 


(76) 
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(77) 


^1  M.(c)V 


-4,0') 

n,  m 


a 


A.U)  A’V) 

n,  m  t  n,m 


n  I?  i*  A,(j)  /  A,(/'k2  A,(/') 

F\dB„(c)F\  dBK(c)Vn,m  =  (<*«,«)  V„,  m 


(78) 


i7i  *  i7i  A,  O')  /  A,(/)x 2  A,(/') 

Fl  d«*(c)^laB,(c)<Pn,m  =  (<*„.*)  <Pn>*  • 


(79) 


Solutions  of  these  systems  may  be  computed  numerically,  with  some  effort,  or  determined  through 
the  computation  of  a  collection  of  multidimensional  integrals.  These  integrals  stem  from  the  coordinate 
mapping  'T  from  the  array  -4  and  the  segment  of  the  unit  sphere  S 2  on  which  we  know  the  scattered  field 
E ;  i.e.,  the  push  forward  ¥ :  A  — »  S2  is  a  diffeomorphism  so  that  the  pullback  T  -1  exists  and  'T  _1 :  S2  —> 
A.  In  this  sense,  the  measured  signal  on  the  two-dimensional  array  admits  the  representation 


oo  n 

£=  s  z  (so* 

n  =  1  m  =  —n 

The  next  result  tells  us  what  the  analysis  basis  for  the  measured  electric  scattering  amplitude  is  on  A; 
namely,  and  V  also  provides  a  more  evident  formulation  of  the  eigenvalue  problem, 

which  must  be  solved  to  determine  the  singular  values  }  and  . 

Proposition  3.  Let  A  be  a  smooth  finite  two-dimensional  manifold  with  a  boundary,  such  that  'T  :  A  -A  S2 
and  'T_1  :  S2  — >  A  are  the  usual  push  forward  and  pull  back  mappings.  Moreover,  let  E^  be  the  measured 
electric  scattering  amplitude  on  the  array  A  and  let  F |0B  (c)  and  F \  dBK(c)  denote  the  restriction  of  F  to 
dBR(c)  and  its  Hilbert  adjoint.  Then 


4,(1) 

VI/ 

t  n,  m 


=  v^yf0'? 


-1 


A, (2)  -  ^ 


(81) 


(82) 


and  a 


4,0) 

n,  m 


J  = 


1 , 2  may  be  determined  from 
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(83) 


=  (<i,))2vSIC'rl 


(84) 


Proof.  Let  E *  denote  the  electric  scattering  amplitude  on  the  unit  sphere  S2  and  the  array  A,  respectively, 
and  let  j  denote  the  magnetic  current  giving  rise  to  the  two  fields.  Then, 


L  (A) 


(85) 


L\A) 


(86) 


=  o',Fr«,(c)^i:#'p'l>iW 


(87) 


/•  (1)  4,(1\ 

0»  Gn,  m*Pn,  m  \2(R*) 


(88) 


=  (PJ, 


a 


4,(1)  4,(1) 
n,m  Tn,  m 


V(/?V 


(89) 


Since  P jj  is  the  projection  of  any  vector  field  in  L2(R\  C 3) 
’ this  amounts  to  saying 


onto  the  orthonormal  basis 


rV  4,(1) 

F\  dBK(c)V„m 


=  a 


4,(1)  4,(1) 
t  n,  m 


(90) 


holds  for  Y| /^’,(„l) 


V^,>  '  .  Similarly,  examining  the  coefficients  K 


in  the  much  the  same  way. 


K,m  = 


(A) 


(91) 
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(92) 


=  <^l 

=  (/.FIWaV^1)^  (93) 


-  /;  4,(2)  4,(2). 

k/>  ® n,  m  ^Pn,  m  ) 


i-V) 


(94) 


=  <P2y. 


A,  (2)  4,(2). 

®n,m  *P n,m  'z,2(/?3)’ 


(95) 


implies 


F\  *BR{c)\f 


-4,(2) 

n,  m 


a 


A,(2)  A,(\) 
n,  m  Tn,m  9 


(96) 


and,  more  precisely,  that  \|/*’^2)  =  x  a  V52P|"ovP  1  .  This  completes  the  proof. 

The  main  result  presented  in  this  section  owes  itself  in  part  to  what  is  known  as  Picard’s  Theorem. 
The  theorem  essentially  provides  a  representation  for  a  compact  linear  operator  A  between  two  Hilbert 
spaces  H\  and  H2  in  terms  of  the  operator’s  singular  system,  as  well  as  a  means  to  assess  whether  a  given 
element  of  the  second  space  H 2  is  also  an  element  of  the  closure  of  the  range  of  A.  We  take  a  moment  to 
state  the  theorem  and  refer  to  [8]  for  its  proof. 

Theorem  3  (Picard).  Let  A  :  H\  — »  //2  be  a  compact  linear  operator  from  the  Hilbert  space  H\  into  the 
Hilbert  space  H2  with  the  denumerable  singular  system  {0„,  <|)n,  \|/n}”=  , ;  i.e.. 


(97) 


and 


A  Vn  = 


(98) 


and  let  (.,.)  denote  the  inner  product  on  H2.  Then,  the  equation  Af  =  g  is  solvable  if  and  only  if 

*  1 

g  e  N(A  )  and 
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oo 


<  oo  . 


(99) 


I 


CT 


2 

n 


Moreover,  any /of  the  form 


f=  I 

n  =  1 


(g,  v„) 


4>n- 


(100) 


solves  Af=  g. 

Finally,  we  state  the  following  range  characterizations  of  F,  which  we  will  need  for  the  proof  of  the 
made  theorem  to  follow. 

Lemma  1.  Let  SR  (c X2)  be  of  two  spheres  of  radii  Rj  and  R2  having  centers  Cj  and  c2,  such  that 
Br  (cx)  contains  BR  (c2) .  Then 


(101) 


2 

Proof.  Let  BR  /cj  2)  be  as  stated  above.  Let  fx  2  e  TL  ( dBR  c^x  2))  •  Then  the  trivial  extension  of  the 
form 


j  f\(x),  xedBR(Ci ) 

{  0,  xedBRi(c{) 


(102) 


proves  the  claim. 

Next,  we  demonstrate  that  in  the  case  when  the  closure  of  BR  (c,)  has  no  intersection  with  the 
closure  of  BR  (c2) ,  then  the  only  common  element  in  the  ranges  of  F|3gs  (c  )  and  F|5g  ((. )  is  the  trivial 
electric  scattering  amplitude  E M  =  0.  This  amounts  to  saying  that  two  disjoint  conducting  scatterers  can 
never  produce  the  same  electric  scattering  amplitude  on  any  open  set  away  from  their  support. 


Proposition  4.  Let  BR  (c,  2)  be  as  before.  Then 


^/?,(ci)  n  Br2(c2^  ~  ^ 


R(F | 


=  {0}- 


(103) 
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Proof.  Each  tangential  component  of  E on  S2  satisfies  Rellich’s  lemma;  c.f.  [8].  Hence,  if  we  treat  each 
component,  this  is  a  immediate  consequence  of  corollary  5.1.2  in  [9], 

Proof  of  the  Main  Theorem.  Application  of  Picard’s  Theorem  to  the  identified  spaces 
H\  =  TL  ( dBR(c ))  and  H2  =  TL  (S  )  and  the  range  characterizations  of  F|3B  (r)  now  prove  the 
theorem  in  the  previous  section.  Specifically,  Green’s  theorem  proves  1)  =>  2),  while* the  asymptotics  for 
the  generalized  Fourier  coefficients  for  any  free  solution  of  the  Maxwell  system  on  the  complement  of 
Br(c)  proves  2)  =>  3).  Lastly,  assuming  3)  and  knowing  that  the  range  of  F |dB  (t;)  is  dense  in  TL2 (S2)  and 
employing  lemma  1  and  proposition  4  we  arrive  at  1).  (C.f.  [9]  for  the  details  of  these  arguments.)  This 
completes  the  proof. 

Given  Picard’s  Theorem,  and  our  range  characterization  of  the  operator  F|5b  ( (  ),  we  now  have  a 
definitive  test  that  can  determine  whether  a  scatterer  is  fully  within  some  region  El  of  interest  by  means  of 
testing  the  convergence  of  the  sum 


ini2 


11 7T2(aQ) 


z 


n  =  1 


n 


z 


m  =  -n 


<£~>hO 

2 

<C(Q> 

(104) 


In  theory,  if  the  sum  does  not  converge,  then  we  can  conclude  that  the  scatterer  is  not  fully  within  the 
test  region  El,  whereas  if  the  sum  does  converge,  we  can  conclude  that  the  scattering  support  of  the  scatter 
is  fully  contained  within  this  set. 


We  note  that  this  formulation  states  both  Theorems  1  and  2  at  the  same  time;  namely,  £2  may  play  the 
role  of  the  sphere  BR(c )  or  any  other  finite  multiply  connected,  possibly  nonconvex  set.  In  the  more  general 
setting,  once  we  have  determined  the  appropriate  singular  system  for  the  operator  F  restricted  to  the 
boundary  of  £2 — i.e.,  we  solve  the  eigenvalue  eigenfunction  equations  in  proposition  3  with  BR(c)  simply 
replaced  by  the  more  general  El — we  may  perform  the  same  Picard  test  and  conclude  whether  the  scatterer 
is  fully  within  El  or  not,  as  indicated  by  the  finiteness  of  the  series.  There  is  good  reason  to  believe  that  this 
test  may  be  done  in  the  electromagnetic  case  discussed  here,  again  on  the  basis  of  the  numerical  procedure 
and  results  obtained  in  [7]. 
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5.  SUMMARY  AND  CONCLUSIONS 


The  results  presented  here  are  fundamental  to  many  advanced  radar  detection  and  imaging  efforts 
and  offer  us  an  avenue  by  which  the  three-dimensional  shape  of  rotating  bodies  may  be  determined  with 
narrowband  monostatic  sensors. 

We  have  rigorously  developed  and  presented  the  underlying  theory  of  the  scattering  of  a  single 
monochromatic,  or  narrowband,  incident  electromagnetic  wave  and  its  interaction  with  a  perfectly 
conducting  object.  Moreover,  we  have  carefully  detailed  how  this  scattering  problem  leads  to  a  viable  and 
robust  reconstruction,  or  three-dimensional  imaging  algorithm,  that  is  capable  of  estimating  the  size  and 
shape  of  either  a  convex  or  multiply  connected  nonconvex  scatterer  with  this  extremely  limited 
information. 
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